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Abstract. 

Stringent physical constraints relate the J/ip produced at a given transverse momen- 
tum pt to the Bjorken-cc of the initial gluons. We present a new approach which takes them 
into account in order to explicitely investigate the p^-dependence of the shadowing effect on 
the J I ifj production. Using the J/ tp rapidity and pr spectra extracted from = 200 GeV 
p + p data from PHENIX, we build a Glauber Monte-Carlo code which includes shadowing 
in two alternative ways: multiple scattering corrections or Q 2 evolution of parton densities. 
We present our results in d + Au collisions at the same energy, notably providing the first 
prediction of the J /ip nuclear modification factor as a function of pr, and compare them to 
the available data by adding some nuclear absorption effect. 
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1. INTRODUCTION 

Relativistic heavy ion collisions may be used as a tool to produce the quark-gluon 
plasma (QGP), the state of the hadronic matter for extreme temperature and 
density [1]. The J/i/j(cc) is foreseen to be a sensitive probe to the QGP formation, 
due to effects such as dissociation by colour screening in the deconfined medium [2] . 
The PHENIX experiment at RHIC recently measured the J/ip production in 
Au + Au collisions at ^s^n = 200 GeV [3]. The yield is quite suppressed with 
respect to p + p collisions [4] scaled by the equivalent number of binary collisions. 
But the interpretation relies on a proper subtraction of the cold nuclear matter 
(CNM) effects, known to impact on the J/ip production in an ordinary nuclear 
environment - such as in any p(d) + A collisions where the QGP can not be 
created. J/ production in d + Au collisions at ^/snn = 200 GeV was measured by 
PHENIX [5] to establish the experimental baseline. At RHIC, we shall focus on two 
CNM effects: (i) the shadowing (initial-state effect) due to the modification of the 
free nucleon structure function by the nuclear environment, and (ii) the breakup of 
correlated cc pairs (final-state effect) due to multiple scattering with the remaining 
nucleons from the incident nuclei, often referred to as "nuclear absorption" . 

Intensive theoretical work has been carried out about nuclear effects on the 
structure functions [6], especially on the shadowing (see the recent review [7]) and 
its effects on the cc production (see e. g. [8, 9, 10]). In this contribution, we present 
our approach to explicitely investigate the dependence of the shadowing on the 
J/ip transverse momentum p?. Two different shadowing models (CF and EKS) 
will be used to get the J/ip production cross-section in nuclei collisions; the latter 



being described in the framework of a Glauber [11] Monte-Carlo. We will extend 
the prior work reported in [12] by adding the nuclear absorption in order to be able 
to compare the obtained CNM effects to PHENIX d + Au data, notably providing 
a very first prediction of their pr-dependence. 



2. THE ^-DEPENDENCE OF SHADOWING MODELS 
2.1. Shadowing observables 

The free nucleon structure function FJ? is the sum over the various parton species 
of their momentum distributions xfi(x,Q 2 ) weighted by their charge e, squared: 
F^ = J2i e "i ■ x fi{ x -iQ 2 ) where i stands for each of the parton species, i.e. all valence 
and sea (anti)quarks, and the gluons, Bjorken-x is the fraction of the nucleon 
momentum carried by the parton, Q 2 is the energy scale of the process used to 
probe F2 i and /j is the parton density (PDF). 

Various processes can be used: for e. g. electroweak processes like the deep 
inelastic scattering (DIS) [13] I + N — > I + X, or the Drell-Yan (DY) process 
p + N —> fi + fi~X, both sensitive to the quarks and antiquarks PDFs. They provide 
indirect contraints on the gluon PDF - notably via the deviations of F<^ from 
the Bjorken scaling caused by gluon radiation. At RHIC, J/ip production mainly 
proceeds through gluon fusion [14], hence probing the gluon PDF. 

The bound nucleon structure function is studied in I + A or p(d) + A thanks to 
the same processes. A summary of the available measurements in the (x,Q 2 ) plane 
used to determine the nuclear PDFs (nPDF) can be found in [15]. At a given Q 2 , 
the ratio Rp 2 (x,Q 2 ) of the bound to the free nucleon structure functions deviates 
from unity: the shadowing corresponds to the small-x region (usually x < 0.1) 
where Rp < 1, while the ant i- shadowing region with Rp 2 > 1 lies at intermediate-x 
(usually 0.1 < x < 0.3). So any process with initial nuclear partons in these regions 
will be suppressed (resp. enhanced). If shadowing is the sole effect on the J/ip 
production, then only the gluonic part Rf of Rp 2 is needed to build the correction 
factor that relates the cross-section in p(d) + A to the one measured in p + p: 



a 



V{A)A = Shadow X(AU>^ (1) 



where -R^ adow = f(R^) and (Y coll ) is the average number of collisions in p(d) + A. 



Eq. (1) is identical for p + A and d + A since shadowing effects in the deuterium 
are negligible. Only its spatial extension must be taken into account. 

As can be intuited from Eq. (1), the relevant experimental observable is the J/ijj 
nuclear modification factor RdA, defined as the ratio of the production cross-section 
in d + A to the one in p + p scaled by (N co \\), or equivalently computed as: 



(N coll )dN^/dy 



u-iv dA /u,y 



where dN^^ pp y/dy are the measured J/ip yield per rapidity unit. 



2.2. Physical origin and models of shadowing 



Shadowing appears as the consequence of coherence effects [6, 7]. 

CF approach. Let us describe the shadowing in the target nucleus rest frame. 
In DIS, the incoming virtual photon fluctuates into a qq pair long before reaching 
the nucleus. At high energy (or at low-x), its coherence length 1 lc = l/2m^x can 
become of the order of the nuclear radius, leading to a coherent interaction with 
several nucleons at once. The cross-section per nucleon is then reduced, which gives 
birth to shadowing. The hadronic component of the virtual photon interacts with 
a pomeron "emitted" by a nucleon, and shadowing is the outcome of the Gribov 
theory [7, 8, 16] composed of a multiple-scattering (multi-pomeron exchanges) 
formalism and a diffraction component from interactions between the pomerons 
(which essentially are gluons). The CF shadowing model [8, 9, 12] belongs to such 
approaches. Following the Schwimmer unitarization scheme [17] to sum diagrams 
with triple pomeron interaction, the correction factor due to shadowing in p + A 
collisions can be expressed at fixed impact parameter 2 b as: 

def 1 

Ri^(b,y, PT ) = RL(b,y, PT ) = l + ATA[h)F{y ^ T) ( 3 ) 

where y is the center of mass rapidity of the produced particle, pt is its transverse 
momentum, Ta(6) is the nuclear profile function related to the Woods-Saxon 
distribution z) at z longitudinal location in space by f dzpA{b,z) = AT^b). 
The function F(y,px) accounts for initial interactions between gluons. It is given 
by the integral of the ratio of the triple pomeron cross-section to the single 
pomeron one: 



y max x d 2 a PPP 



rVmax 

F(v,Pt) = 4tt / dy— ■ 

J Vmin 



dydt 



C [exp (A . y max ) - exp (A . y min )] (4) 



t=o 



where t = {p — p') 2 is the usual variable in DIS related to the difference between 
the incoming and outgoing quadri-momentum of the nucleon, y min = In ( ^^p- J 

and y ma x = f hi TV with y > (y < 0) for the projectile (target) hemisphere, 
with Ra being the nuclear radius, s the square of the center-of-mass energy per 
collision and the transverse mass of the produced meson, = y'm^ +Pt- 
C is a function of the parameter A, the pomeron-proton coupling <7™(0) and the 

triple pomeron coupling r ppp (0), both evaluated at t — 0: C = gpp ^ A — — . The 
values of C and A can be fixed from data on DIS. We have used C = 0.31 fm 2 and 
A = 0.13 as in [8] and references therein. 



mjv stands for the nucleon mass. 

b denotes a position in the transverse plane. 



EKS approach. Equivalently, the physical picture of the (anti-)shadowing can 
be viewed from the Breit frame (nucleus infinite momentum frame) . Low-x partons 
spread over a large longitudinal distance Az, proportional to l/(rriNx). Partons 
from different nucleons may spatially overlap, interact and fuse, hence increasing 
the parton density at higher-re at the expense of the one at low-x (conservation of 
the nucleon momentum). The EKS model makes use of the fact that nPDFs and 
PDFs are different as a starting point, without addressing the physical origin. This 
kind of models rather use the experimental data to provide a parametrization 3 of 
all the parton-specie dependent ratios Rf(x,Ql) of nPDF/PDF, at some fixed scale 
Ql, which enter into Rp 2 (x,Ql) by a sum. Although perturbative QCD (pQCD) 
cannot give the absolute Rf(x,Q 2 ), it can predict their evolution with Q 2 starting 
from some initial Qo > Aqcd ; thanks to the DGLAP [18] equations. 
The EKS model is named after the EKS98 parametrization [19] used to compute 
the shadowing correction factor. In EKS98, the ratios Rf(x,Q 2 ) are obtained at 
Q 2 = 2.25 GeV 2 , are evolved at LO for Q 2 < 10 4 GeV 2 and are valid for x > 10" 6 . In 
this model, the spatially-dependent shadowing in p + A collisions is given by [10]: 

< ad ow(M,Q 2 ) = Rf(b^Q 2 ) % i + N p j^^ [R?(x,Q 2 ) - 1] (5) 

where pA(b,z) is normalized to get J d 2 b J dzpA(b,z) = A, N p is a normalization 
factor chosen to have (l/A) J d 2 bj dz pA(b,z)Rf(b,x,Q 2 ) = Rf(x,Q 2 ). There is 
a simple idea behind Eq. (5). At large b, the nucleon density is small, so we 
expect the target nucleons to behave as free nucleons: the shadowing effects 
should be negligible. At 6 = 0, we rather probe the center of the target nucleus, 
where the vicinity of a high density of nucleons should lead to larger effects. 
In this equation, the integral over z includes the target nucleus material that 
the projectile nucleon traveled through, by its longitudinal path at an impact 
parameter b. So the spatially-dependent shadowing is obtained by assuming that 
the projectile parton interacts coherently with all the target partons localized 
within a cylinder, its axis being defined by its longitudinal path and its transverse 
section area by the nucleon transverse area a^.. The average nucleon density being 
Po = (l^rr 3 )" 1 = 0.17 nucleon/fm 3 , we get = nr 2 = 3.94 fm 2 . Therefore, the 
number of target partons which contribute to shadowing will be larger at small b. 



3 However, deriving the nPDFs and hence the ratios Rf(x,Q 2 ) from the data is difficult: in the 
nuclear case, there is an additional dependence in A and Z, and there is no DIS data below 
x < 5.1CT 3 at Q 2 > Aq CD . So the nuclear ratio R^(x,Q 2 ) for the gluons is less constrained at 
low-x, which is accessible by high energy colliders for heavy quark production. Note also that 
the existing DIS data do not allow to determine the ^-dependence of the ratios Rf(x,Q 2 ). But 
we know from PHENIX results that there is a centrality dependence of i?dAu, the J/ip nuclear 
modification factor. 



2.3. Where is the px dependence? 



For the J/ip, the shadowing is nothing more than a production process where 
the initial partons are picked within the nPDFs, the latter exhibiting a different 
dependence on x and Q 2 with respect to the PDFs. The underlying production 
process puts stringent physical constraints on the initial partons that can make 
a J/ip at some given rapidity y and transverse momentum p?, through quadri- 
momentum conservation and partonic cross-section dependence on s = snnX\X2- 
Hence, specifying the production process comes down to defining the shadowing 
dependence on y and px- The simplest production process is 

g + g-^cc (6) 

At Pt = 0, quadri-momentum conservation results in: 

xi,2 = ^=exp(±y) (7) 

V S AW 

Keeping this simplest process for pt ^ implies that the initial partons cannot be 
colinear: they have to carry an intrinsic transverse momentum, later on transferred 
to the created cc. If the px of the J/i/j has such an intrinsic origin, then: 

x\,2 = —j^=ex-V\±y) (8) 

The available theoretical predictions of the EKS shadowing [10] were all made 
at pt = with the use Eq. (7). In our work [12], we will investigate how these 
predictions, made at RHIC energy, are affected by the introduction of a non-zero 
Pt when using Eq. (8). We also add the non-zero pt in the energy scale: 

Q 2 = (2m c ) 2 + (p T ) 2 (9) 

where m c = 1.2 GeV/c 2 is the charm mass value, in accordance with [10]. 

About the CF shadowing [8], the explicit p T -dependence of Eq. (3) rely on the 
assumption in Eq. (8). The previously published predictions [9] are made at a fixed 
value of y and a unique value of pt ~ since they only consider niT = 3.1 GeV. In 
our work [12], we will allow both y and Pt to vary within their full spectra. 



3. J I if) PRODUCTION IN OUR MONTE-CARLO 

To investigate the CNM effects on the J/i/j production, we implemented a Monte- 
Carlo framework, with three main steps. 

Step 1: nucleus-nucleus collisions described with a Glauber [11] model. For each 
A + B collision (an event), the value of the impact parameter b is randomly chosen 
relative to the so-called "minimum bias" distribution (2irbdb). Then, for each 
nucleus, the nucleon positions are randomly chosen according to the nuclear density 



profile, defined by the Woods-Saxon [20] parametrization for any nucleus A>2 and 
by the Hulthen [21] parametrization for the deuterium. Within this code, we can 
determine the number of target nucleons in the path of each projectile nucleon and 
calculate an event-by-event N co \\. Such a binary collision is considered to occur if 
the distance d between two nucleons satifies red 2 < <r NN , where a^N stands for the 
nucleon-nucleon cross-section (at RHIC, ct^n = 42 mb at ^s NN = 200 GeV). 

Step 2: kinematics for the produced J ftp candidates. For each nucleon-nucleon 
collision, a J/ip candidate can be produced (with an arbitrary production cross- 
section aj/i/j < ^-rp-), with y and px randomly chosen in the respective input 
spectra. The latter are given by fits to the recent PHENIX p + p data [4] at 
^/s NN = 200 GeV. The angular orientation if of px in the {p X) p y ) plane is also 
random and uniformly distributed in [0,27r]. The Bjorken-x carried by the initial 
partons are then computed with Eq. (8), and the scale Q 2 with Eq. (9). In order 
to remain in the physical phase space domain, we require that < xi,X2 < 1. 

Step 3: involving CNM effects. To account for shadowing, only some J/i/j 
candidates are randomly allowed to become "real" J/ip, since the production cross- 
section of the real J/tp is aj/^ corrected by the factor R^ dow , where 

- rL shadow - rL shadow -""shadow V lu / 

and ^sha^i are gi yen by Eq. (3) and (5) for the CF and EKS models, respectively. 
To account for the nuclear absorption, any Monte-Carlo J/ip that is "wounded" 
by the remaining incident nucleons is tagged as broken. The cross-section <7brcak-up 
is used to define the distance d required to get a J/V-nucleon interaction: red 2 < 
(Tbreak-up- Finally, the J/ip nuclear modification factor Rab is computed as followed: 

dN no t w ounded real J/i/)/dy , , 

R AB= JT- T J-, I 11 ) 

aMj/il> candidate/ U-V 



4. RESULTS AT RHIC ENERGY AND DISCUSSION 

Fig. 1 shows -RdAu as a function of y when considering or not a ^-dependence of 
the CF and EKS models. Here, shadowing is the sole effect studied, so -RdAu can be 
identified to -R^ adow - Adding an intrinsinc p^-dependence has a slight effect only 
(more visible for the EKS case), because of the average value of pt in Eq. (8): 
(pr) < 2 GeV/c < mj/^ at RHIC energy. For CF, F{y,pT) is a monotonic function, 
increasing with y and decreasing with px- So at fixed y, a larger pt results in 
a smaller value of F(y,pT) and hence a larger R^ adow from Eq. (3). For EKS, a 
larger pt leads to a larger Q 2 . In the higher-x part of the anti-shadowing region, 
this leads to a smaller Rf(x,Q 2 ) > 1, and hence a smaller -R^ adow from Eq. (5) 
(the contrary is obtained for the lower-x part of the shadowing region). Fig. 1 also 
shows that the amount of antishadowing is larger for EKS compared to CF. 

Fig. 2 presents the expected -RdAu as a function of (iV co ii) for the three rapidity 
windows accessible by the PHENIX detector. There is a remarkable difference 



between the effect predicted by CF and EKS models at backward rapidity: the 
smaller amount of anti-shadowing for the CF model shows up here. This does 
not come as a surprise, since the CF model was originally built to describe the 
coherence effect that leads to the depletion of the nuclear structure function at 
low-a;. Its validity domain for the present energies starts around y > —2. 




FIGURE 1. J/ip nuclear modification factor in d + Au collisions as a function of rapidity for 
CF (top) and EKS (bottom) shadowing models. Three input px distributions are used: pt = 0, 
p T from \y\ < 0.35 and pr from 1.2 < \y\ < 2.2. 
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FIGURE 2. J/ip nuclear modification factor in d + Au collisions as a function of the number 
of collisions for CF (left) and EKS (right) shadowing models for three rapidity windows: backward 
(—2.2 < y < —1.2), central (|y| < 0.35) and forward (1.2 < y < 2.2) regions (from up to down). 

Our main result lies in Fig. 3, which shows the expected i?dAu as a function of px 
for both models. The pr-dependence is rather significant: it can lead to amplitude 



variations as large as about 20%. The EKS model exhibits a stronger dependence 
with Pt- At backward rapidity, EKS et CF notably show opposite behaviours. 
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FIGURE 3. J/tp nuclear modification factor in d + Au collisions as a function of the J/ip 
transverse momentum for CF (left) and EKS (right) shadowing for the backward, central and 
forward rapidity regions (from up to down). 




FIGURE 4. J I ih nuclear modification factor in d + Au collisions as a function of rapidity 
for CF (top) and EKS (bottom) shadowing on top of the nuclear absorption and compared 
to PHENIX data [5]. Bars stand for point to point uncorrelated errors (both statistical and 
systematical uncertainties) and boxes for the point to point correlated systematical errors. 

The predictions with shadowing as the sole effect overshoot the data, which 
advocates for an additional nuclear absorption effect. It is parametrized by the 
break-up cross-section Ubreak-up, not calculable from first QCD principles. For each 



shadowing model, we choose 0break-up in order to get the best description of -RdAu 
vs y as measured by PHENIX [5] . The results are shown on Fig. 4: different values 
of the cross-section are required for the CF model (0" brcak _ up = 1.5 mb) and the 
EKS model (obreak-up = 2.8 mb i. e. the same value as in [5]). The corresponding 
expectation of -RdAu as a function of pt is compared to PHENIX data [5] on 
Fig. 5. Unfortunately, their large uncertainties and limited range in px do not 
allow to draw firm conclusions. At backward rapidity, their trend with seems 
to disagree with the EKS model. For the other rapidity regions, both models give 
an approximate description of the data. At forward rapidity, they both seem to 
underestimate the slope, which may indicate some missing effect (like a possible 
broadening of the transverse momemtum of the initial gluons due to the multiple- 
scattering experienced as the incident nucleon travels through the nucleus). 
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FIGURE 5. J/ip nuclear modification factor in d + Au collisions as a function of the J/ip 
transverse momentum for CF (left) and EKS (right) shadowing on top of the nuclear absorption 
and compared to PHENIX data [5] for the three rapidity regions. Bars stand for point to point 
uncorrelated errors and boxes for the point to point correlated errors. 



5. CONCLUSIONS AND OUTLOOK 



In summary, we compared the J/ip shadowing expected from two models, CF and 
EKS, in d + Au collisions at RHIC described in the framework of a Glauber 
Monte-Carlo. In general, the suppresion is stronger for CF, which implies a smaller 
break-up cross-section for the nuclear absorption effect needed to match the data. 
We investigated their p T -dependence based on the assumption that the J/ift's Pt 
would be of intrinsic origin and we used the y and pt spectra from the p + p data 
as an input of the Monte-Carlo. Although the resulting -RdAu vs y or vs (N co \\) 
was only slightly affected when introducing such dependence, this approach lead 
to the first predictions of -RdAu as a function of pr- Due to large uncertainties and a 
limited range in p T , the available data do not allow to state definitive conclusions. 



However, the data seem to have an increasing trend with px at backward rapidity, 
which sound to be in accordance with the expectations from CF rather than EKS. 
Hopefully, there is much higher statistics d + Au data recently taken at RHIC. 

For the models, there are two major sources of uncertainty: the parametrization 
of the ratios nPDFs/PDFs used in EKS-type models (see [22, 23]), and the way 
the J /if) acquires its pr- The production process g + g — * cc + g, with strictly coli- 
near initial gluons, can also give a non-zero pt (extrinsic origin). The corresponding 
pr-differential cross-section from [24] reasonably matches both RHIC and Tevatron 
data. The next step will be to use this computation to relate J / ip's (y,pr) to (x±, x 2 ) 
in the extrinsic scheme. The main advantage will be to free the Monte-Carlo from 
any input spectra from the p + p data in order to get predictions at LHC energies. 
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